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Abstract. We report the discovery of a hitherto undetected population of Wolf-Rayet stars in the young galactic 
open cluster Westerlund 1. Optical spectroscopy of the cluster identified 11 such objects; provisional classihcation 
suggests that 6 are nitrogen rich (WN) and 5 carbon rich (WC). Including the previously identified Blue, Yellow 
& Red Super- & Hypergiants, Westerlund 1 clearly has a very rich population of massive post-Main Sequence 
objects. To date, the post-MS population of Westerlund 1 is significantly larger than that of any other galactic 
young open cluster - with the possible exception of the Arches - implying that it is potentially amongst the most 
massive young clusters yet identified in the Local Group. 
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1. Introduction 


2. Observations 


The highly reddened young open cluster Westerlund 1 
( hence forth Wd 1) was first identified by Westerlund 
( 1961 ). Subsequ ent b roadband ph otome tric surveys by 
Borgman et al. ( jl 970 ), Lockwood (1974) and Koornneef 
( 1977 ) suggested the presence of a number of both early 
and late type supergiants, while a comprehensive photo¬ 
metric and spectroscopic survey of the b right est cluster 
members was presented by Westerlund ( 1987 ; West87). 
Despite reporting the presence of a large number of very 
luminous (L >10 5 Lq ) transitional objects only one fur¬ 
ther (photometric) study of the cluster has been made 
(Piatti et al. 1998). 

Recently, radio continuum observations of Wd 1 re¬ 
vealed that a number of the cluster members appeared to 
be as sociated with very bright radio sources (Clark et al. 
199S, Dougherty, Clark & Waters, in prep.). Motivated 


by these results we obtained low resolution optical spec¬ 
troscopy of a number of the brighter cluster members in or¬ 
der to provide an accurate spectral classihcation for them. 
In this paper we present the first results of this program; 
the discovery of a significant population of Wolf-Rayet 
(WR) stars within the cluster. 


Send offprint requests to : J. S. Clark, e-mail: 
j scOstar.ucl.ac.uk 

* Based on observations collected at the European Southern 
Observatory, La Silla, Chile (ESO 67.D-0211) 


Spectroscopy of cluster members over the red/near- 
IR spectral region (~6000-11000A) was taken on 2001 
June 23-25 from the ESO 1.52-m telescope at La Silla 
Observatory, Chile. The telescope was equiped with the 
Loral #38 camera and the #1 (night 1) and #13 (night 
2 and 3) gratings, giving dispersions of ~ 5 A/pixel and 
~ 7 A/pixel - leading to resolutions of ~ llA and « 16A - 
respectively. Data reduction was accomplished with pack¬ 
ages within the Starlink software suite. 

Due to the crowded nature of the held, each long slit 
integration included a number of different cluster mem¬ 
bers. Examination of the fainter objects present in several 
of the exposures revealed the presence of a number of ob¬ 
jects with rich emission line spectra (see Figs. 1 and 2). 
Given that the integrations were optimised to avoid satu¬ 
rating on the brighter cluster members, the serendipitous 
sources are of a low S/N ratio - though sufficient to iden¬ 
tify the emission line objects as a previously unidentihed 
population of massive, hydrogen depleted WRs. 


3. Results 

Despite the low S/N of many of the spectra, it is imme¬ 
diately possible to identify both nitrogen rich WN (6 ob¬ 
jects) and carbon rich WC (5 objects) stars (Figs. 1 and 
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2 respectively); a finding chart and co-ordinates for each 
object are presented in Fig. 3 and Table 1. 

Accurate determination of the spectral types of the 
WN and WC stars in this spectral region is difficult, given 
that most commonly used diagnostics lie at shorter wave¬ 
lengths. However initial spectral classification of th e WR 
cand idates using the catalogues of Vreux et al. ( 1983 , 


1990J) was possible. 


For preliminary classification of the WC candidates 
we use the ratio of the Ciii(8500 A)/Civ(8856 A) 
and Cn(9900 A)/Ciii(9710 A) lines (Vreux et al. 1983| ; 


Howarth & Schmutz 1992; Crowther priv. comm.). 


Civ(8856 A) is absent from the spectra of all candidates 
with an upper limit to emission in candidates F, E and 
C constraining spectral types to later than WC6. No con¬ 
straints are possible for H and K due to weak or no emis¬ 
sion in Cm(8500 A). Cn(9900 A)/Ciii(9710 A) ratios of 
0.17±0.04, 0.18±0.03 and 0.18±0.04 for candidates F, E 
and H respectively, suggest WC9 classifications while a 
ratio of <0.08 suggests a WC8 classification for candidate 
Cfl while a classification for candidate K is not possible. 

Classification of the WN spectra proved more difficult, 
as there are fewer lines present and there is no linear pro¬ 
gression in e.g. the N lV:He II ratios that might be used for 
classification (Vreux et al. 1983| ). Based on the strength 
of the N iv 7103-7128A feature we can exclude extreme 
WN3/WN9 subtypes for all objects, since these are the 
only spectral types for which it is not seen in emission. 
Candidate A appears to be a WNE (WN4-5) given the 
large lines widths (~3000 kms^ * 1 ) and lack of strong Hei 
7065A emis s ion (s een for all stars between WN6-8; Vreux 
et al. 1983 , 1990| ). On the basis of emission in this line 
and the small emission line widths, candidates G, I and D 
appear to be WN6-8 objects. Classification of candidates 
B and J is complicated by the low S/N in the blue regions 
of the spectra due to the CCD response curve, but are also 
probably WNL objects. 


4. Discusion & Conclusions 


Unfortunately, uncertainty in the distance and redden¬ 
ing estimates to Wd 1 (West87; Piatti et al. 1998 ), cou¬ 
pled with the poorly determined bolometric corrections for 
many massive transitional objects make determination of 
the luminosities of the WRs and the other evolved cluster 
members difficult. Indeed, the lack of accurate luminos¬ 
ity estimates for the yellow hypergiant candidates (YHG; 
West87) is particularly concerning, given that bolometric 
luminosity is one of several classification criteria for such 


objects (e.g., de Jager et al. 1998). If we are to constrain 
the post-MS population of Wd 1 in order to determine 
such fundamental properties as cluster age and mass we 
must address these issues. 


1 EW(CII)/EW(CIII)=0.05-0.08 for WC8 (from WR 135, 
113 and 53) and EW(CII)/EW(CIII)=0.14-0.3 for WC9 (from 
WR 92 and 104) 



Fig. 1 . Spectra of the newly discovered WN candidates in 
Wd 1, with prominent transitions identified. 



Fig. 2. Spectra of the newly discovered WC candidates in 
Wd 1 with prominent transitions identified. 


For the WR candidates A, E and F (for which West87 
provide photometry), adopting the reddening and distance 
estimates given by West87pj and the (conservative) bolo- 


2 note that radio observations of the WR stars suggest that 
the cluster is more distant than the alternative distance of 

1 kpc proposed by Piatti et al. 
prep. 


1998 


see Dougherty et al., in 
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Fig. 3. R broadband finding chart for the newly discovered Wolf Rayet candidates in Wd 1. Note that the exact 
counterparts of candidates H, J & K should be regarded as provisional since the crowded field in these regions of the 
cluster make identification of the correct counterpart difficult. Candidates A, E & F correspond to stars number 72, 
241 & 239 respectively in the notation of West87. 


Table 1. Co-ordinates (J2000) for the newly identified 
Wolf Rayet stars in Wd 1 determined from 3.6cm radio 
images (Dougherty, Clark & Waters; in prep.) Formal er¬ 
rors are o a = ±0.003 s , and <j$ = ±0.04", although given 
that the crowded fields in the vicinities of candidates H, 
J & K make identification of the correct counterpart dif¬ 
ficult, the errors for these objects are likely to exceed the 
formal error quoted. 


Candidate 

a 

S 

A 

16h47m8.324s 

—45°50'45.51" 

B 

16h47m5.354s 

—45°51'05.03" 

C 

16h47m4.395s 

—45°51'03.79" 

D 

16h47m6.243s 

—45°51'26.48" 

E 

16h47m6.056s 

-45°52'08.26" 

F 

16h47m5.213s 

-45°52'24.97" 

G 

16h47m4.015s 

-45°51'25.15" 

H 

16h47m3.905s 

—45°51'19.88" 

I 

16h47ml.668s 

—45°51'20.40" 

J 

16h47m0.885s 

-45°51'20.85" 

K 

16h47m2.697s 

—45°50'57.35" 


luminosities comfortably exceed 10 5 Lq. For luminous A, 
F and G stars the strength of the OI 7774A feature can 
be used to provide an additional measure of the absolute 
visual magnitude^. Comparison of the YHG candidates 
W4, W12, W16 and W265 and the YSG W7 (notation 
from West87) to the luminosity dine strength calibrations 
of Ferro & Mendoza ( 1993| ) and Slowik & Peterson ( 1995[ ) 
indicates that they are all likely to be intrinsically highly 
luminous. Indeed the EW(Ol) for these objects are sig¬ 
nificantly greater than any of the stars in these studies - 
including the the bona fide YHG p Cas (L*=10 5 ' 7 Lq) - 
implying luminosities of L* 10 5 L S . Additionally, high 
mass loss rates for these objects are implied by radio detec¬ 
tions (W4, W12, W16 and W265) and broad Ha emission 
(W7, W12, W16 and W265); therefore all these objects 
meet the classification criteria of de Jager (1998) for bona 
fide YHGs. 


metric corrections of Smith et al. (1994) we find that their 


3 Note that our values of EW(Ol) are consistent with those 
given by West87; Negueruela & Clark (in prep). 


















4 


J. S. Clark and I. Negueruela: Wolf-Rayet stars in Wd 1 


4.1. Comparison to other clusters 

Including the WRs and YHG candidates, the large popu¬ 
lation of luminous post-MS objects in Wd 1 (West87 and 
Table 4) suggests that it is unique in both the number 
and variety of massive post-MS objects present. Of the 26 
clusters within the solar circle studied by Eggenberger et 
al. ( [2002 ), 6 contain both B and RSGs. Of these, only 3 - 
Collinder 228, Trumpler 27 and Berkley 87 - also contain 
WRs (1-2 per cluster; Table 4) and only Trumpler 27 con¬ 
tains a yellow SG, albeit of significantly lower lumin osity 
(~10 4 7 Lq) than those in Wd 1 (Massey et al. 200 1| ; note 
however that they claim the cluster is not co-eval.). 

Of the Galactic Centre clusters, the Quintuplet clus¬ 
ter appears to contain 8 WRs of both WN and WC types 
along with a number of early OB supergiants and a sin¬ 
gle RSG (1O 4 ' 9 L 0 ; Figer et al. 199Sp . In addition to the 
single luminous RSG IRS 7 (10 ^T©; Carr et al. 200CQ , 
Paumard et al. (2001) detect 16 “helium stars” in the 


Galactic Centre cluster, of which they suggest the 7 nar¬ 
row line objects correspond to mass losing BSGs (possible 
LBVs) and the 9 broad li ne obj ects WRs (noting that a 
further 3 stars Blum et al. ( 1996 ) suggest are WCs lie out¬ 
side their f-o-v). Only the Arches cluster - which is signifi¬ 
cantly younger - appears to contain a comp arable number 
of WRs to Wd 1, with Blum et al. ( 2001 ) identifying 15 
candidate 04 If/WN7 stars on the basis of narrow band 
imaging. However, as with the Quintuplet and Galactic 
Centre clusters, the rich population of very luminous cool 
stars present in Wd 1 is absent [jj 

Clearly, such comparisons indicate that Wd 1 is both 
very young and very massive; however uncertainties in the 
post-MS evolution of massive stars, exacerbated by the 
lack of an identifiable MS turnoff and accurate bolometric 
luminosities for the evolved stars makes determination of 
the age and total mass of Wd 1 difficult. 

4.2. The age and mass of Wd 1 

Follow ing th e analysis of the Quintuplet cluster by Figer 
et al. ( 1999| ), the presence of WC stars - apparently the 
most evolved stars present in Wd 1 - implies a lower limit 
to the age of 2.5 Myr. The presence of a number of very lu¬ 
minous RSGs within Wd 1 potentiall y pro vides an upper 
limit to the cluster age; Figer et al. ( 1999 ) suggests that 
the Quintuplet cluster requires an age of >4 Myr given the 
presence of a single (low luminosity) RS G, br oadly con¬ 
sistent with the age (7 Myr) Carr et al. (200C) claim for 
IRS 7. However, large uncertainties in the mass loss rate 
for very luminous cool stars render estimates of their ages, 
lifetimes and progenitor masses highly uncertain; partic¬ 
ularly concerning given the large population of YHGs 
within Wd 1. 


4 the very young clusters NGC 3603 and R 136 also contain 

a few high luminosity WRs although these are thought to be 

very massive O stars where the high mass loss rates simulate 

the spectra of more chemically evolved lower mass WRs (e.g. 

de Koter et al. 1997). 


Considering their extreme rarity, the lifetime of YHGs 
is probably less than 10 5 * * * yr for any luminosity and pro¬ 
genitor mass. At very high luminosities it is likely that a 
very large mass loss rate limits the YHG to a single pas¬ 
sage from red to blue across the HR diagram, resulting in 
a short lifetime (<30,000 yr; Stothers & Chin 1999Q . At 
lower luminosities dynamical instabilites in the outer at¬ 
mosphere of the star result in multiple blue loops for the 
star out of the RSG region, leading to a longer lifetime as 
a luminous yellow star (e.g. Stothers & Chin 2001). Such 
estimates are qualitatively consistent with the results of 
unpublished simulations for the behaviour of Mj n i t i a i=25 
& 40 Mq stars (Maeder & Nieuwenhuijzen priv. comm. 
2002) which suggest YHG phases after 6.9 & 4.4 Myr last¬ 
ing ~49000 and 2700 yr respectively; i.e. the YHG phase 
occurs earlier and is shorter the more massive the progen¬ 
itor is. 

Despite the many uncertainties, the present stellar 
population of Wd 1 appears consistent with an age of order 
4-8 Myr if the cluster is co-eval, suggesting it is potentially 
younger than previo usly thought (7 and 8±3 Myr; West87 
& Piatti et al. 1998 , respectively). 

The uncertainties in post-MS evolution and incomplete 
stellar census inevitably makes a determination of the to¬ 
tal cluster mass uncertain. Maeder & Meynet (1994) find 
that for solar metallicities, adopting twice the standard 
mass loss rate results in the appearance of a WN phase 
for stars >25 M 0 - consistent with the findings of Massey 
et al. (2001) - and a WC phase at >40 Mq. Given that 
the super- and hyper-giants are less chemically evolved 
than the WRs, their progenitors are likely to have been 
less massive. 

Regarding the completeness of our sample we estimate 
that the spectral survey of potential cluster members of a 
similar visual magnitude as the fainter WR candidates is 
at best 33 per cent complete, while preliminary analysis 
of our low resolution optical spectroscopy (Negueruela & 
Clark, in prep.) finds many additional supergiant candi¬ 
dates, e.g., W70 and W71 (BSG), W32 and W33 (YSG) 
and W75 and W237 (RSG). Additionally, we might expect 
that if originally present very m assive stars will have been 
lost to SN; Maeder & Meynet (|1994| ) suggest that stars 
of >85 Mq will have a lifetime comparable to the lower 
estimate of the cluster age. 

Nevertheless, conservatively assuming that the present 
stellar census is complete (see Table 4) and furthermore 
that the progenitor masses for these objects were >30M 0 
we can derive a lower limit to the initial cluster mass. 
Adopting a Salpeter mass function ( N(M ) oc M~ a ) with 
upper and lower cut offs of 100 M 0 & 0.2 M 0 respectively 
and a slope, a=2.35) we might expect a total mass of stars 
of ^750 Mq for every star with an initial mass of ~30 M 0 , 
leading to a mass estimate of a few xlO 4 M 0 (S. Goodwin 
priv. comm. 2002). 

This lower limit to the initial mass of Wd 1 suggests 
that it is directly comparable to the galactic centre clusters 
such as the Arches cluster (4xl0 4 M 0 ; Portegies Zwart et 
al. 2001). More reasonable estimates of completeness and 
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Table 2. Transitions and equivalent widths for the 
WN candidates. Errors are estimated at 10% for lines 
longwards of ^7000A and 20% for those shortwards. 
Transitions for which the S/N is too poor to attempt an 
identification are indicated with ‘S/N’. 


Candidate 

Transition 

A 

G 

I 

B 

J 

D 

Hen 6560 

27 

S/N 

9 

S/N 

S/N 

S/N 

He i 6678+Hen 6683 

- 

30 

37 

S/N 

S/N 

S/N 

Civ 7061+Hei 7065 

- 

59 

50 

45 

- 

46 

Niv 7103-28 

110 

B1 

B1 

B1 

19 

B1 

Hen 8237 

33 

20 

19 

12 

7 

11 

Hell 5-8? 

38 

- 

- 

- 

- 

- 

Hen 10124 

106 

56 

42 

33 

32 

41 


progenitor mass suggest a mass for Wdl of a few x 10 5 Mq 
making it by far the most massive young Galactic cluster, 
and one of the most massive in the Local Group. This 
conclusion is further reinforced by considering the large 
number of YHGs within Wd 1, which is compar able t o 
the total population of the Milky Way (6; de Jager 1998| ). 
Assuming a lifetime for the YHG s of ~ 25000 yr, following 
the arguments of Geballe et al. ( 2000| ) for the likelihood 
of finding several examples of a short lived evolutionary 
phase within a single cluster also leads to the conclusion 
that a very large O star population (^several hundred) is 
required to produce the number of YHGs observed. 

Evidently a combined spectroscopic and photometric 
approach to both identify the MS turn-off and to properly 
classify evolved stars will be required to accurately deter¬ 
mine the age and mass of Wdl. However if the above esti¬ 
mates are correct then Wd 1 would appear to be a Galactic 
equivalent of the super star clusters observed in merging 
and interacting galaxies and may possibly be more mas¬ 
sive than the 30 Doradus cluster in the LMC. Therefore, as 
well as providing a unique laboratory for studying hot star 
evolution, Wd 1 would provide an unprecidented insight 
into an extreme mode of cluster formation, previously not 
thought to be occuring in the Milky Way. 


Table 3. Transitions and equivalent widths for the WC 
candidates; terminology as for Table 2. 


Candidate 

Transition 

F 

E 

C 

H 

IC 

Hen 6560+Cn 6578 

42 

49 

S/N 

S/N 

13 

He i 6678 

9 

17 

S/N 

S/N 

- 

Cn 6725-42+C hi 6727-73 

B1 

B1 

S/N 

S/N 

- 

Cn 6780 

28 

45 

S/N 

S/N 

- 

Cm 7037 

B1 

B1 

S/N 

S/N 

- 

He i 7065+Cn 7064 

15 

32 

S/N 

S/N 

- 

Cn 7236+C hi 7210-12 

55 

108 

45 

24 

3 

Cm 7772-96 

S/N 

8 

S/N 

- 

- 

Cm 8196 

12 

15 

16 

5 

- 

Cm 8256 

B1 

9 

B1 

- 

- 

Cm 8328-59 

39 

27 

18 

10 

- 

Cm 8500 

14 

17 

11 

6 

- 

Cm 8664 

19 

12 

- 

- 

- 

Hen 9225+ Cn 9224 

16 

25 

20 

14 

- 

Cm 9705 

103 

139 

97 

50 

26 

Cn 9903 

18 

25 

S/N 

9 

- 


Table 4. Numbers of B & RSGs and WRs for galactic 
open clusters containing one or more of each type of ob¬ 
ject. The numbers of B & RSGs for Wd 1 are those from 
West87; preliminary analysis of our data suggest they 
are likely to be lower limits. Note that the Quintuplet 
cluster also contains 5 ‘Cocoon’ objects with feature¬ 
less IR spectra that have been proposed as dusty WCL 


stars. (“Massey, DeGioia-Eastwood & Waterhouse |2001 
and references therein; b Figer, McLean & Morris M. 1999 


“Paumard et al. |2001 
|2000t e Blum et al. fcOOll ). 


Carr, Sellgren & Balachandran 


Cluster 

Age(Myr) 

Ybsg 

NrsG 

Nwr 

Wd 1 

4-8 

>6 

>2 

>11 

Trumpler 27“ 

NA 

8 

1 

2 

Collinder 228“ 


3 

1 

1 

Berkley 87“ 


4 

1 

1 

Quintuplet 6 

4±1 

14 (+ 2 LBVc) 

1 

8 

Gal. Center 6 ’ d 

3-8 

7 He i 

1 

9 He i 

Arches 6 

2-4.5 

0 

0 

15 
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